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Abstract 

The common endosymbiotic Wolbachia bacteria influence arthropod hosts in multiple ways. They are mostly recognized for 
their manipulations of host reproduction, yet, more recent studies demonstrate that Wolbachia also impact host behavior, 
metabolic pathways and immunity. Besides their biological and evolutionary roles, Wolbachia are new potential biological 
control agents for pest and vector management. Importantly, Wolbachia-based control strategies require controlled 
symbiont transfer between host species and predictable outcomes of novel U/o/bacWa-host associations. Theoretically, this 
artificial horizontal transfer could inflict genetic changes within transferred Wolbachia populations. This could be facilitated 
through de novo mutations in the novel recipient host or changes of haplotype frequencies of polymorphic Wolbachia 
populations when transferred from donor to recipient hosts. Here we show that Wolbachia resident in the European cherry 
fruit fly, Rhagoletis cerasi, exhibit ancestral and cryptic sequence polymorphism in three symbiont genes, which are exposed 
upon microinjection into the new hosts Drosophila simulans and Ceratitis capitata. Our analyses of Wolbachia in 
microinjected D. simulans over 150 generations after microinjection uncovered infections with multiple Wolbachia strains in 
trans-infected lines that had previously been typed as single infections. This confirms the persistence of low-titer Wolbachia 
strains in microinjection experiments that had previously escaped standard detection techniques. Our study demonstrates 
that infections by multiple Wolbachia strains can shift in prevalence after artificial host transfer driven by either stochastic or 
selective processes. Trans-infection of Wolbachia can claim fitness costs in new hosts and we speculate that these costs may 
have driven the shifts of Wolbachia strains that we saw in our model system. 
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Introduction 

The ubiquitous intracellular a-proteobacterium Wolbachia 
pipientis manipulates reproductive biology of many arthropod 
species in order to warrant its own transmission in host 
populations (reviewed by [1,2]). Wolbachia are maternally inherited 
and favor infected females by inducing reproductive phenotypes 
with cytoplasmic incompatibility (CI) as the most common ([3-5]). 
Besides reproductive host manipulations, Wolbachia can also affect 
nutritional and metabolic pathways of hosts ([6,7]), host develop- 
ment and lifespan (reviewed in [8-1 1]), provide protection of hosts 
from pathogens and parasites ([12-15]), as well as affect host 
mating behavior and facilitate host speciation ([16,17]). 

Wolbachia have attracted major attention as potential novel 
biological control agent against the increasing threat that insect 
populations pose to human health and agriculture either by 
vectoring pathogens, or by directly causing damage to crops and 
life stock ([18]). Wolbachia could suppress insect populations 
through the induction of CI in a way analogous to the Sterile 
Insect Technique (SIT; reviewed in [19]), manipulate vectorial 



capacity through host lifespan-shortening effects ([11]) or vector 
refractoriness to pathogens in mosquitoes ([10,15,20-23]). Wolba- 
chia could also be used in combination with the widely used SIT 
([24]) that has encountered some problems with male fitness, 
mating performance, sperm size and number due to gamma 
irradiation of individuals ([25,26]). However, all of these Wolbachia 
applications require that (t) Wolbachia strains are transferable 
between diverse host systems and insect genera, (a) the transferred 
strains are capable of inducing the expected phenotype such as CI, 
life shortening or resistance against pathogens, and (uij transferred 
infection and their expected phenotypes persist stably in the novel 
host. So far, several authors have reported successful transfer of 
Wolbachia by microinjection from donor to recipient hosts within 
the same insect order, followed by confirmation of their phenotype 
in the novel host ([27-29,11]). The third requirement of 
phenotypic, and thus genomic stability has not yet been tested 
extensively, although artificially transferred symbiont strains can 
potentially experience phenotypic changes ([30]). Genotypic 
changes might include point mutations and genomic rearrange- 
ments triggered and facilitated by symbiont infection dynamics 
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upon arrival and successful establishment in the novel host. In 
addition, pre-existing variability of Wolbachia in the donor host 
such as the presence of spurious genomic polymorphism in 
neotropical Drosophila species ([31]) and tsetse flies ([32]), as well 
as the existence of low-titer multi-infections ([33]) might affect the 
outcome of the artificial transfer and the stability of the expected 
phenotype in the new host. Recent studies have uncovered 
phenotypic plasticity of Wolbachia over evolutionary short time 
periods within the same host, and also upon transfer between 
different host backgrounds. Weeks et al. ([34]) revealed a rapid 
switch of wBA from a parasitic ([35]) towards a more benign 
mutualistic state in Californian populations of Drosophila simulans. 
Adaptation of the symbiont to a novel host cell environment 
resulted in significant phenotypic changes of wMelPop when 
transferred between hosts ([36]). In Drosophila melanogaster wMelPop 
causes early death via over-replication in mainly nervous tissue 
([37]). Trans-infection experiments caused the same lifespan 
reduction in D. simulans and Aedes aegypti ([38,11]). For the latter 
host, aMelPop was initially pre-adapted to Aedes cell lines, before 
its successful establishment in mosquitoes via embryonic trans- 
infection. Re-introduction of the Aedes-cell line adapted Wolbachia 
from cell lines into their native host D. melanogaster, however, 
resulted in fluctuations of maternal transmission efficiency, lower 
titers, and a reduced life shortening effect ([36]). These findings 
implicate that both host and symbiont interact dynamically and 
co-evolve rapidly within relatively short time periods. 

It is so far unknown how genotypic integrity of Wolbachia is 
affected by recombination, genetic drift or selection after artificial 
host transfer. Here, we have monitored Wolbachia genome 
dynamics and population structure experimentally by utilizing 
two different Wolbachia-insect host species, D. simulans ([29]) and 
the Mediterranean fruit fly Ceratitis capitata ([28]). Both hosts were 
microinjected with Wolbachia of the European cherry fruit fly 
Rhagoletis cerasi ([29,28]). We tested if the transfer of Wolbachia, in 
particular of the strain «)Cer2, (i) induced de novo structural 
rearrangements, and/or (it) enhanced sequence polymorphism 
within the symbiont genome post microinjection. The potential 
appearance of additional Wolbachia subvariants, hereby designated 
haplotypes, that are distinctive from the reference sequence of the 
symbiont in its native host ([39]) can either be assigned as de novo 
mutations or ancestral haplotypes that had already persisted in the 
donor at low frequencies and thus had escaped earlier detection. 
We hypothesized that the structure of the original Wolbachia strain 
used for trans-infection plays a crucial role in the performance 
within the new host. It is frequently assumed that Wolbachia strains 
are monoclonal populations, and thus, only de novo mutations will 
contribute to potential diversification of the symbiont in the 
recipient host. In contrast, an alternative situation of Wolbachia 
strains that represent diverse populations of haplotypes such as in 
native tsetse flies ([32]) at varying frequencies would allow for 
detection of sequence variation after host transfer that was not 
triggered by new mutations. The trans-infection event may solely 
shift haplotype frequencies and thus enable detection of rare 
ancestral haplotypes. 

To assess the potential for new structural rearrangements of 
wCer2 in its two new host species, we took advantage of the 
unusual high numbers of mobile genetic elements in Wolbachia 
genomes ([40-45]) with their capacity to trigger insertions, 
inversions as well as ectopic recombination ([46-50]). Our 
analysis, however, did not reveal rearrangements. To test for 
sequence polymorphism, we analyzed SNP (single nucleotide 
polymorphism) frequency of three Multi Locus Sequence Typing 
(MLST) Wolbachia genes and traced a>Cer2 sequence heterogeneity 
in original and new hosts. In the course of this in-depth analysis we 



found in microinjected D. simulans clear signs of an unexpected 
infection with wCerl of the original host that had remained 
undetected for over 150 generations. We then discuss whether 
heterogeneity in new hosts is caused by ancestral Wolbachia 
sequence polymorphism or arises through new mutations. 

Materials and Methods 

2.1 Insect lines 

Fly stocks of D. simulans and the Mediterranean fruit fly C. 
capitata known to be infected by wCer2 were used in this study. 
They had been microinjected with cytoplasm from the Wolbachia- 
infected cherry fruit fly R. cerasi thirteen ([29]) and eleven years ago 
([28]). Rhagoletis cerasi is naturally multi-infected with up to five 
strains, roCerl - wCtrb ([51,39]). Based on sequence analysis of 
MLST genes, «;Cerl, tyCer2 and <x>Cer4 are A supergroup strains, 
wCcrb a B supergroup strain, and a>Cer3 a recombinant strain 
([39,52]; http://pubmlst.org/wolbachia/). For the first host- 
transfer experiment, embryos of Wolbachia-Ccee D. simulans 
(Noumea TC, generated by tetracycline treatment over three 
generations; [53]) were injected with egg cytoplasm from Austrian 
R. cerasi donors in 2000 ([29]). From this experiment, six GO 
isofemales resulted in wCer2 infected lines RC20, RC21, RC33, 
RC45, RC50 and RC78 that had to be further selected for 
Wolbachia in consecutive generations. For selection, DNA was 
extracted from multiple females and infection status was 
determined via Wolbachia-speciiic PCR. Offspring of females that 
tested positive for Wolbachia was used to proceed into the next 
generation via sibling mating ([29]). This selection regime was 
paused between G20 (2001) and G140. In 2007, selection for 
Wolbachia was continued after only five out of the six initial 
isofemale lines, RC20, RC21, RC33, RC45, and RG50 tested 
positive for «;Cer2. The Wolbachia strain ojCerl, however, was not 
detected since it was considered lost between Gl and G2 ([29]). 
For the second host-transfer experiment, embryos of the C. capitata 
Benakeion line were injected with Wolbachia from R. cerasi from 
Austria and Italy (Sicily) in 2002 ([28]). This resulted in two 
infected C. capitata lines ([28]) and one of these, WolMed88.6 
harboring 0)Cer2, was included in our study. DNA extracts of 
approximately G50 post-infection were kindly provided by K. 
Bourtzis' laboratory (University of Ioannina, Greece). Wolbachia- 
free D. simulans Noumea TC and D. melanogaster w ( Wof l '" s ) were 
used as negative controls, Wolbachia-infected DSR (D. simulans 
Riverside, California; [54]) and the D. melanogaster Harwich strain 
(WoP os ) were used as positive control. All Drosophila lines were kept 
on standard medium at 24°C. 

Antibiotic treatment of donor and recipient hosts. For 
Wolbachia depletion, embryos of R. cerasi were transferred to larval 
medium containing a final concentration of 0.02, 0.025 and 0.05% 
(w/v) tetracycline and incubated at 24°C until reaching the third 
instar (L3). Antibiotic larval media and R. cerasi individuals were 
kindly provided by K. Koppler from the Center for Agricultural 
Technology Augustenberg (LTZ), Stuttgart, Germany. D. simulans 
RC20 and RC50 were placed on standard Drosophila diet 
containing 0.03% tetracycline for two consecutive generations 
before the presence of Wolbachia was tested. 

2.2 Molecular isolation and characterization of Wolbachia 
from R. cerasi and their novel hosts 

DNA extraction. For PCR, cloning and sequencing, high 
quality genomic DNA was extracted from individual pupae or 
adults using the Puregene DNA Purification Kit. For Southern 
blots, genomic DNA was extracted from individual adult flies and 
processed following the protocol from [55]. For Wolbachia 
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depletion assays, DNA was extracted from tetracycline-treated 
adults of trans-infected D. simulans RC lines using the Puregene 
DNA Purification Kit (Qiagen). DNA was stored at -20°C until 
use. 

MLST-PCR, cloning, sequencing and sequence 
analysis. We analyzed the frequencies of SNPs in the three 
Wolbachia MLST genes gatB (WD.0146, aMel), coxA (ITO0301, 
oMel) and ft$Z (WD_0723, aMel) in donor and recipient hosts 
([56,57]). General coxA and JisZ primer sets were used as in [56]. 
For gatB, an additional primer set amplifying a 404-bp fragment, 
was designed (gatBY 5'-gatttaaatcgtgcaggggtt-3' and gai£_450R 5'- 
ttgaattaaatcaattttatcctgg-3'). To selectively target o)Cerl, we used 
a strain-specific primer set described in [39] plus the VNTR-141 
primer set from [48]. For all PCR reactions a Biometra T3000 
Thermocycler was used. PCR products were purified using the 
peqGOLD Gel Extraction Kit, inserted into the pTZ57R/T 
vector and then transformed into competent DH5a Escherichia coli 
cells. Clones containing the insert were Sanger sequenced at the 
University of Chicago Cancer Research Center (UCCRC-DSF). 
Sequences were analyzed using ApE plasmid editor vl.10.4 (M. 
Wayne Davis), CodonCode Aligner Version 2.0.3 (CodonCode 
Corporation) and the BLAST algorithm (www.ncbi.nlm.nih.gov). 

For D. .simulans recipients, we analyzed SNP frequencies for each 
line separately (RC20, RC21, RC33, RC45, and RC50) plus for 
the pool of all trans-infected lines (RC) in order to maximize 
sample size. The automated base-calling in CodonCode Aligner 
software detected SNPs in many clones from single individuals. In 
order to verify their authenticity, we visually inspected all SNPs in 
the corresponding chromatograms from both forward and reverse 
reads. All ambiguous SNPs were excluded from the final data set. 
Confirmed SNPs were then divided into two groups: recurrent 
SNPs and true singletons. Recurrent SNPs refer to nucleotide 
positions that were detected in independent clones in either 
different lines of new hosts, in R. cerasi clones only, or in both 
systems. Unique SNPs that were found in single clones only, but 
appeared reliable in the sequence chromatogram, were classified 
as singletons. 

Anticipating that PCR accuracy was strongly impacted by the 
performance of the enzyme polymerase, we first determined the 
error base line of the Taq DNA polymerase used for all assays. We 
did not use a proof-reading enzyme as the Promega Go Taq® 
Flexi DNA Polymerase used in all experiments is one of the best 
non-proof-reading high quality and high performance polymerases 
on the market (see Promega notes available at www.promega. 
com). Based on re-PCR and re-sequencing of known coxA, fts£, 
and gatB fragments inserted into the pTZ57R/T cloning vector 
from independent batches of this polymerase (data not shown), we 
calculated the following error base lines: 0 SNPs in 4.44 kb oicoxA 
(10x444 bp), 0 SNPs in 2.90 kb of/fe? (10x290 bp), and 2 SNPs 
in 6.00 kb of gatB (14x429 bp). Compared to the published mean 
estimate for standard non-proof-reading Taq DNA polymerases of 
0.21 SNPs/kb ([58]), our Promega Go Taq® Flexi DNA 
Polymerase control assay thus resulted in lower DNA polymerase 
error rates (0.15 SNPs/kb). We deposited 33 coxA,ftsZ, and gatB 
sequences, which represent rare allelic variants and nonsense 
mutations at GenBank (accession numbers KF17541-17573). 

Restriction Fragment Length Polymorphism (RFLP) 
analysis via single fly Southern hybridization. We deter- 
mined the structural integrity of the bacterial chromosome in the 
novel hosts via RFLP-analysis with highly dynamic Wolbachia 
marker sets: Insertion Sequence elements (IS) and Variable 
Number of Tandem Repeats (VNTRs). 0.5-1 |lg genomic DNA 
from single flies was digested with 6U Hindlll (New England 
Biolabs, USA) for 4 hrs at 37°C. After high-resolution vertical gel 



electrophoresis ([55]), gels were vacuum-blotted onto a positively 
charged nylon membrane (Hybond™-N+, GE Healthcare, UK). 
Membranes were hybridized with [a- 32 P]dCTP-labeled specific 
probes of IS and VNTR loci. Probes were prepared with the 
RediprimeTM II DNA labelling kit (GE Healthcare, UK) and 
exposed to high sensitivity X-ray films (Kodak, Germany). Probe 
primers were designed with respect to the annotated genome of 
aMel of D. melanogaster (NCJ302978; [40]). For transposon 
probing, three repeats greater than 200 bp, belonging to different 
IS families (IS3, IS5, and ISNew) in the aMel genome of D. 
melanogaster were chosen ([40]). VNTR probes were VNTR-141, 
consisting of tandemly repeated 141-bp units, located at coordi- 
nates 89,003-90,332 in aMel ([48]) and VNTR-144, consisting of 
11.8 copies of a 144 bp repeat unit located at 34,727-37,210 in 
aMel (MR, unpublished). /fedlll-digested DNA of the lambda 
phage (New England Biolabs, USA) was utilized as size marker. 
Size of fragments was determined with respect to this standard, 
allowing the comparison of number and size of fragments between 
donor R. cerasi and new hosts D. simulans and C. capitata. 

2.3 Ovary Screening Assay 

We estimated the general fecundity cost of the artificially 
generated infection by comparing the ovaries of recipient and 
uninfected control lines of D. simulans. Analysis was performed 
between G168 and G182 post microinjection and followed the 
fecundity assay by [59]. Fertilized females were raised on standard 
food and were dissected in sterile PBS ten days after eclosion. 40 
ovaries per recipient Dmsophila line were screened. Fecundity status 
of ovaries was estimated according to number of mature eggs in 
the ovary: 0 eggs = class I; 1-2 eggs = class II; 3-9 eggs = class III; 
1 0 and more eggs = class IV. Only eggs at stage 1 4 of oogenesis 
([60]; indicated by dorsal appendages) were counted. 

2.4 Statistical analysis 

For statistical analysis SPSS 16.0 and GraphPad Software 
(www.graphpad.com) were used. SNP frequencies were analyzed 
using % with Yates Correction (2x2 contingency table); two tailed 
P-values indicated significant differences between values when 
<0.05. To detect potential traits of positive or negative selection, 
synonymous to non-synonymous substitutions per site were 
calculated using SNAP (Synonymous and Non-synonymous 
Analysis Program) provided at http://www.hiv.lanl.gov. This 
program uses the Nei-Gojobori corrected path counting method 
that adjusts for counts via Jukes-Cantor plus the weighting of 
pathways from one codon to another according to an equi- 
probable model for each possible codon-to-codon path ([61]). 

Results 

3.1 Conserved Wolbachia genome synteny after artificial 
transfer into novel recipient hosts 

We analyzed the genome synteny of artificially transferred 
o,'Cer2 infection from R. cerasi into D. simulans and C. capitata with 
five marker probes for RFLP mapping (Figure S2). Although IS 
and VNTRs had earlier been reported as hypervariable entities of 
Wolbachia genomes, we did not detect any structural re-arrange- 
ments for the five tested loci in roCer2 of the novel hosts (Figure 
S2 and Data SI). 

3.2 SNP frequency in gatB, coxA, and fcZof recipient host 
populations 

Prior to all sequencing experiments, we determined the base line 
error rate during polymerase chain reaction of the Taq polymerase 
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used in our lab. Consequently, we analyzed SNP frequencies in 
three Wolbachia loci (gatB, coxA, and ftsZ) in trans-infected D. 
simulans and C. capitata. For gatB, we sequenced and analyzed 
29.1 kb: 12x404 bp from each of RC20, RG21, RC33, RC45, 
RC50 of D. simulans, plus 12x404 bp from WoMed88.6 of C. 
capitata (Table SI). In total, 38 SNPs randomly dispersed along 
the 404 bp gatB gene fragment (Table 1) were detected. Out of 
these 38 SNPs, 6 (16%) were found recurrently, and 32 (84%) were 
singletons. Recurrent SNPs were SNP- 11 (2x); SNP-42 (4x); 
SNP-93 (2x); SNP-186 (2x); SNP-250 (2x), and SNP 253 (2x 
ochre mutation, i.e., CAA to TAA, occurring in RC of D. simulans 
and in WoMed88.6; Table 1). Overall SNP frequency for gatB of 
»Cer2 was 1.01 SNPs/kb (Table 2A). For coxA of roCer2 we 
analyzed 16.4 kb (37x444 bp amplicon) with an average SNP- 
frequency of 1.64 SNPs/kb (Table 2B). For fts^, we analyzed 
7.6 kb of «;Cer2 (16x478 bp amplicon) with an average SNP- 
frequency of 0.78 SNPs/kb (Table 2C). Overall SNP-frequencies 
for all three genes were significantly higher, i.e., 5-1 lfold, than the 
error base line we had determined before for the Tag polymerase 
(compare frequencies of 0.78, 1.01, and 1.64 SNPs/kb to 0.15 
SNPs/kb from Taq polymerase). 

3.3 Uncovering cryptic co-infection with the wCerl strain 

Unexpectedly, we detected sequence traces of roCerl in three 
lines of microinjected D. simulans, although this strain had 
previously been considered as lost from the microinjected lines 
between Gl and G2 ([29]; and Figure 1A). To confirm the 
presence of wCerl independentiy, we performed PGR analyses 
utilizing primer sets specifically targeting wsp of either oderl or 
wCer2, respectively ([39]; Figure IB). According to the signal 
intensity of the wsp amplicon during electrophoresis, RC20 
harbored «;Cerl at high densities whereas «)Cer2 could no longer 
be tracked (Figure IB). Wsp sequence reads from RC20 revealed 
exclusively the a)Cerl haplotype (Figure 1C). In RC33 and 
RC45, we traced «jCerl co-infecting wCer2 (Figure IB). RC50 
showed no signals of wCerl but of a single infection with o)Cer2 
(Figure IB). In addition to wsp, the presence of wCerl in RC20 
and its absence in RC50 was confirmed by the diagnostic VNTR- 
141 locus via PGR (Figure ID). 

After verification of the wCerl sequences we analyzed the SNP 
frequency of this strain using gatB and coxA in a small sample of 
new hosts. We tested 1,212 bp for gatB with an average SNP- 
frequency of 2.48 SNPs/kb (Table 2A). For coxA we sequenced 
4.9 kb (3x305 bp amplicons from RC20 and 9x444 bp from 
!4WMed88.6 and RC20 and RC33; Table 2B) with an average 
SNP frequency of 1.63 SNPs/kb. SNP frequencies for JtsZ of 
MjCerl in novel hosts were not determined due to the low coverage 
(Table 2C). 

Regarding the heterogeneity detected in wCtr, we hypothesized 
that the SNPs can (i) either represent ancestral, hidden sequence 
polymorphism, i.e., multiple Wolbachia haplotypes already present 
in the donor host, or (iij have arisen de novo following 
microinjection. To test this hypothesis, we compared SNP 
frequencies in the novel hosts with the frequencies in the donor 
of the Wolbachia strains. Since original donor specimens for the 
microinjection experiments into D. simulans and C. capitata were not 
kept as voucher material, we sequenced gatB, coxA, and jtsJ? from a 
broader representative range of other R. cerasi populations across 
Europe instead. 

3.4 SNP frequency baseline in gatB, coxA, and ftsZ of 
donor and comparison with recipient host populations 

We analyzed Wolbachia SNP frequencies from R. cerasi samples 
collected all over Europe ([39]). In total, we sequenced 8.9 kb 



(22 x404 bp amplicons) of «>Cer2 gatB from R. cerasi derived from 
individuals from more than ten collection sites across Europe 
(Table S2 and Table S4). As shown in Table 2A the average 
SNP frequency for gatB of wCer2 from the donor R. cerasi was 
1.01/kb (9 variable positions in 8.9 kb). Comparing these data 
with those of the novel hosts, we did not observe an increase in 
SNP-frequencies per kb of wCer2 gatB in the recipients (Table 2). 
Only one (SNP-42) out of the 38 SNPs, however, determined in 
wCer2 gatB of the novel hosts, occurred also in the original donor 
R. cerasi. 

Since we unexpectedly detected wCerl in three trans-infected 
host lines we included this Wolbachia strain in our SNP analyses 
(however at a lower coverage). For gatB of wCerl, we sequenced 
9.7 kb consisting of a 24x404 bp amplicons data set and 
calculated an average SNP-frequency for gatB of wCerl in R. 
cerasi that was not statistically significant from wCer2 in R. cerasi 
(0.52 vs. 1.01; P= 0.3346; Table 2A). The average SNP-frequency 
in gatB of «;Cerl in novel hosts lines was also not higher than in the 
native host R. cerasi (2.48 vs. 0.52; P= 0.0703), but the small sample 
set (1,212 bp) impeded statistical testing and thus did not allow a 
reliable comparison between both wCer 1 and <x>Cer2. 

For the coxA locus, we sequenced 5.7 kb of oderl (13x444 bp 
amplicons), and 888 bp of roCer2. Average SNP-frequencies for 
i*Cerl and wCtr2 infections were rather low (0.69 SNPs/kb and 0 
SNPs/kb, respectively; Table 2B). It must be taken into account 
that the data set for wCer2 consisted only of two clone reads and 
can thus not be regarded as highly representative. 

For the ftsZ locus, we sequenced 10.5 kb of wCerl (23x478 bp 
amplicons) and 10 kb of wCev2 (21 x478 bp amplicons). At this 
locus SNP frequencies of the two Wolbachia strains were quite 
similar with 1.14 SNPs/kb and 1.10 SNPs/kb for roCerl and 
wCer2, respectively (Table 2C). One SNP in ftsZ of roCerl 
resulted in a transversion from guanine to thymine in the first 
position of the consensus triplet GGA (Gly), thus introducing the 
stop codon TGA to the sequence. We detected this stop codon 
twice and independently in Jts^oi wCevl. 

3.5 Presence of stop codons in n/Cer1 and wCer2 of 
essential Wolbachia genes in both donor and recipient 
host populations 

The recurrent finding of SNPs causing nonsense mutations in 
essential Wolbachia housekeeping genes was highly unexpected and 
hence considered with extreme caution. However, similar to the 
ochre mutation in gatB of wCer2 in D. simulans and C. capitata (see 
3.2), we found additional nonsense mutations in coxA and fts£ of 
wCerl Wolbachia (Table 3). A recurrent SNP (4 x) in ftsZ of twCerl 
of its native host R. cerasi caused a transition of guanine to thymine 
in the first position of a GGA triplet in wCerl (opal mutation). 
Finally we also uncovered an ochre mutation in coxA of wCer 1 in 
recipient line RC20 but as a singleton only (Table 3). In order to 
test whether such Wolbachia pseudogenes might stem from 
translocations onto the host chromosome, tetracycline-treated 
individuals of RC20, RC50, and R. cerasi, plus their corresponding 
untreated controls were tested via gatB PCR for the presence of 
potential nuclear Wolbachia copies. As shown in Figure SI, both 
recipient lines lost the gatB PCR signal after two generations of 
antibiotic treatment, which makes a lateral gene transfer event 
unlikely. Hence alternative scenarios will be necessary for 
explaining these counter intuitive findings, i.e., the persistence of 
nonsense mutations in essential Wolbachia genes (see discussion). 
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wCerl wCer2 wCer4 




Figure 1. Cryptic co-infection with wCerl in wCer2 carrying trans-infected lines. (A) Switch of strain prevalence from wCer2 to wCerl in 
RC20. Asterisk represents time point of line establishment via microinjection in 2000. Generations are indicated on x-axis. (B) Presence of n/Cer1 in 
random samples of RC lines determined via strain-specific wsp PCR. First RC20 sample seems to carry Wolbachia at extremely low density below 
detection limit of wsp PCR whereas the other one gives a bright band with wCerl -specific wsp primer set. DNA extracted from adult Rhagoletis cerasi 
served as positive controls (n/Cer1 and n/Cer1+2). (C) Random 36-bp fragment of the general wsp amplicon showing diagnostic wCer1/wCer2 sites. 
(D) Differentiation between wCerl, n/Cer2 and wCer4 Wolbachia using VNTR-141 PCR. Abbreviations: CC C. capitata, M DNA size marker, T1-3 trans- 
infected RC line sample. 
doi:1 0.1 371 /journal.pone.0082402.g001 



Table 3. Stop codons in gatB, coxA, and ftsZ of wCerl and 
wCer2 Wolbachia. 



Line 


Gene 


Position 


Mutation 


Wo/-infection 








CAA = consensus 


wCer2 


RC20 


gatB 


226/404 


TAA (ochre) 


wCer2 


RC50 


gatB 


253/404 


TAA (ochre) 


wCer2 


Wo/Med88.6 


gatB 


253/404 


TAA (ochre) 


wCer2 








AAA = consensus 


wCerl 


RC20 


coxA 


22/444 


TAA (ochre) 


wCerl 








GGA = consensus 


wCerl 


F37 eastern Sicily 


ftsZ 


25/478 


TGA (opal) 


wCerl 


F38 eastern Sicily 


ftsZ 


25/478 


TGA (opal) 


wCerl 


F40 western Sicily 


ftsZ 


25/478 


TGA (opal) 


wCerl 


F42 western Sicily 


ftsZ 


25/478 


TGA (opal) 


wCerl 



Lane three lists the position of the mutation corresponding to the size of the 
amplified MLST-gene fragment. Lines F37 to F42 represent R. cerasi individuals 
from different populations sampled in Sicily, Italy. 
doi:1 0.1 371 /journal.pone.0082402.t003 



3.6 Trans-infection with wCer is costly for the recipient 

Between G168 and G182 post-micro injection we performed an 
ovary screening assay ([59]) to estimate the costs of establishing an 
artificial infection with Wolbachia. Based on the number of mature 
eggs per ovary, individuals were sorted into the four fecundity 
classes. In total, we screened 40 individuals per D. simulans line 
including 40 individuals from Wolbachia-umnkcted D. simulans 
Noumea TC, which had been used as recipient for the 
microinjection. As shown in Figure 2, percentage of individuals 
in both extreme classes, I and IV, varies highly among new hosts. 
Compared to 80% in the control (Noumea TC), RC21 shows only 
32.5% of individuals in class IV, suggesting poor fecundity 
{P- 0.0004***). Except for RC20, also the other new hosts 
exhibited reduced number of mature eggs (ranging from 52.5 to 
60% in class IV) pointing towards a fitness cost effect related to the 
new symbiont infection. Interestingly, in RC20 percentage of class 
IV ovaries is higher than in all other new hosts and even slightly 
higher then in the control (85% vs. 80%). This might represent a 
positive correlation of fecundity with Wolbachia infection in this 
special case, but would require further testing for verification (see 
discussion). 
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Figure 2. Ovary screen in trans-infected RC lines. Size of ovaries was compared on basis of number of mature eggs in one ovary. Bars represent 
ovary size per line determined for each ovary class: white is class I with no mature eggs; light grey class II (1-2 eggs); dark grey class III (3-9 eggs); and 
black is class IV with 10 or more eggs. Y-axis shows percent of ovaries per class; x-axis shows RC lines plus the Wo/foacWa-unifected D. simulans 
Noumea TC control. Significances based on two-tailed P values from Fisher's exact test are indicated by asterisks. 
doi:1 0.1 371 /journal.pone.0082402.g002 



Discussion 

4.1 No traces of structural re-arrangements detectable by 
means of Wolbachia mobile genetic elements upon 
artificial transfer into novel recipient hosts 

We looked at two highly informative VNTR markers and three 
Wolbachia transposons (IS3, IS5, and ISNew) via Southern 
hybridization, and thus covered at least 30 loci dispersed in the 
wCer2 genome. We did, however, not detect any transpositions in 
these loci. Recent studies demonstrated that Wolbachia carry a high 
percentage of functional transposable elements that can display 
transpositional activity during short term Wolbachia evolution ([47— 
49]). A study of polymorphism in IS element insertions sites and 
VNTRs of wMd revealed that the previously assumed homoge- 
neous Wolbachia infection of D. melanogaster is a set of different 
variants, such as «jMel and wMelCS ([48,62]). However, given 
that so far 20 major families of insertion sequences have been 
classified from 171 bacterial and archeal species ([63,64]) we 
cannot rule out that any transposition effect occurred in our study 
system with any other mobile elements. Thus, as a next step, whole 
genome sequencing of wCer will be the most adequate strategy to 
obtain sufficient information on potential structural rearrange- 
ments in the symbiont genome upon transfer into recipients. 

4.2 Wolbachia strains are most likely bacterial 
populations of diverse haplotypes with shifting 
frequencies 

We determined SNP frequencies for gatB, ftsZ, and coxA, 
between native and recipient hosts and did not detect any SNP 
differences that could have arisen after microinjection. The SNP- 
frequencies in ftsZ, and gatB of wCer2 in the native host R. cerasi 
were with approximately only 1 SNP/kb low. For coxA, the 
frequency was even lower (0 SNPs/kb) but this might not be highly 
representative since we tested a very small sample size. In the 
recipient hosts D. simulans and C. capitata, we determined SNP- 
frequencies in a;Cer2 that were not accelerated compared to R. 



cerasi (see Table 2). Although not statistically significant, we 
observed a trend towards an increase in SNP-frequency, at least in 
gatB of one of the recipients (RC20), at the time when a)Cer2 was 
unambiguously present). We calculated the ratio of the number of 
non-synonymous substitutions per non-synonymous site (djj) to the 
number of synonymous substitutions per synonymous site (ds) for 
gatB of Z£)Cer2, a general indicator for selective pressure acting on 
protein-coding genes. This ratio for gatB, in R. cerasi and RC lines 
(Table S3) corroborated our finding that SNP-frequencies did not 
change upon arrival of Wolbachia in the new host systems. 

Similar to K)Cer2, we found overall low SNP-frequencies in gatB, 
coxA, and ftsZ of iwCerl in the native host R. cerasi, ranging from 
0.5 to 1 . 1 SNP/kb. For the trans-infected hosts, we only analyzed 
the frequency of gatB and coxA, not for JhZ- Both of these data sets 
are rather small and thus still inconclusive. Unexpectedly, we 
revealed cryptic double infections with zeCerl, but this Wolbachia 
strain persisted in lower densities than the predominant <x>Cer2. 
Since the focus of this study was from its onset primarily on wCer2 
and not wCerl, a more extensive analysis of twCerl sequences 
deriving from the trans-infected host systems was not carried out. 
Our coxA sequence data, however, did not indicate an increased 
i*Cerl SNP-frequency in the de novo hosts. In particular, in roCerl- 
gatB, we did not find a statistically significant increase in the SNP 
frequency within RC compared to R. cerasi. 

As shown in Table 1, the variable sites we detected in gatB of 
wCer2 were either present in the native and/or in the trans- 
infected hosts. This finding raises the question whether the 
observed polymorphism represented de novo mutation events or 
ancestral cryptic infection polymorphism. [65] recently showed 
that new Wolbachia haplotypes might be generated by point 
mutations in outer membrane proteins. In our study 68% of the 
SNPs revealed in wCer2 of gatB (26/38) were not detected in the 
donor of the infection but exclusively in the trans-infected hosts, 
suggesting de novo mutation events. Our PCR, cloning and Sanger 
sequencing based approach resulted in relatively low sample size, 
so that this study is not sensitive enough to rule out the existence of 
rare haplotypes in the donor host. Deep-sequencing strategies of 
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donor and recipient with much higher coverage will be of pivotal 
interest to finally uncover potential de novo mutation rates of the 
endosymbiont upon artificial host switch. An alternative reason for 
our failure of detecting shared SNPs in both donor and recipients 
is genetic drift. Through drift effects, allele frequencies change and 
this may result in the loss of certain haplotypes followed by 
consequent reduction of genetic variation within a population. As 
a third explanation, selection for certain, beneficial haplotypes can 
be envisaged. This would lead to subsequent fixation of these 
haplotypes and loss of others. 

We found, however, that SNP-42 in gatB of wCer2 (Table 1) 
occurred recurrendy in both donor and recipient systems (RC20), 
clearly suggesting an ancestral origin. It is possible that the 
mutation in position 42 in the 404-bp gatB fragment represents a 
rare haplotype of wCev2 that coexists with the canonical haplotype 
in the R. cerasi donor ([39]; this study). Moreover, SNP-250 of gatB 
of «)Cer2 that results in the replacement of Arg with Gly (Table 1) 
was found in both heterologous hosts D. simulans and C. capitata, 
but not in the donor. This situation is similar to SNP-253 of gatB of 
wCer2 that results in a nonsense mutation (see below). SNP-11, 
SNP-93 and SNP-186 provide three additional cases for the 
existence of distinctive gatB haplotypes of the i»Cer2 infection since 
all three were repeatedly isolated from hosts that were indepen- 
dently microinjected (Table 1). Hence the twCer2 infection cannot 
be considered monoclonal but a bacterial population of diverse 
haplotypes at varying frequencies. This idea is supported by the 
recent finding from Symula and colleagues who proposed the 
existence of high Wolbachia sequence variation between and within 
individuals of the tsetse fly Glossina fuscipes fuscipes ([32]). 

If wCtv2 infection is a population of haplotypes, rare haplotypes 
within this population might be difficult to detect. Any change in 
the structure of this bacterial population, however, massively 
impacts the frequency of haplotypes. Events that impact the 
population structure as well as population size in such a crucial 
way are referred to as bottleneck events. The artificially transfer of 
wCht2 from its native host R. cerasi into two new hosts, was such a 
bottleneck event and thus manipulated the structure of the original 
wCer2 population. We argue that this resulted in the shift of 
haplotype frequencies in the trans-infected lines. It is likely that the 
polymorphism that we observed in gatB upon arrival in D. simulans 
and C. capitata represented rare haplotypes that already persisted in 
the ancestral and native wCer2 population of R. cerasi and are only 
detectable after the artificially induced bottleneck scenarios. 

4.3 Wolbachia strains accumulate nonsense mutations 
upon arrival in new hosts 

We tested if SNP frequency in Wolbachia genes is increased after 
microinjection, thus suggesting relaxation of purifying constraints 
on these genes. Our results did not explicitly support such an effect 
although a slight trend towards diversifying selection was still 
observed (see 4.2). We revealed, however, that three out of 38 
SNPs (8%) detected in gatB of wCer2 introduced novel pre-mature 
stop codons caused by in-frame ochre mutations, i.e. a transition of 
cytosine to thymine in the first position of a GAA triplet. SNP-226 
was found uniquely in transinfected RC20, whereas SNP-253 was 
found in both recipient hosts independently, i.e. D. simulans (line 
RC50) and C. capitata (line W / o/Med88.6). Novel stop codons were 
not restricted to roCer2 since we also traced them in coxA and JtsZ 
of wCerl. SNP-22 in coxA of roCerl was found uniquely in line 
RC20. In contrast, SNP-45 in JtsZ of roCerl seemed to be of 
ancestral origin, occurring in Sicilian R. cerasi populations only. 

A growing body of empirical evidence has demonstrated that 
Wolbachia genes and even complete genomes are being transferred 
onto insect host chromosomes ([66-69]). Such lateral gene transfer 



events can explain the accumulation of nonsense mutations when 
fully functional copies of these genes are still present in the 
symbiont genome. In order to test for lateral gene transfer, we 
cleared the recipient lines with antibiotics. Symbiont genes 
transferred into the host genome would be not be affected and 
thus still detectable. Our screen of treated RC lines did not 
indicate any gene transfer event, suggesting that the detected 
nonsense mutations are present in cytoplasmic Wolbachia. 

In total, we found two types of nonsense mutations in recipient 
and donor hosts, ochre and opal. In the recipient systems, we 
revealed two in-frame ochre mutations in wCtt2 of gatB, and two 
in-frame opal mutations in roCerl of coxA. In the original donor 
system we determined one in-frame mutation in JtsL of wCerl. 
Mutant tRNA is able of suppressing some stop codons in E. coli 
([70,71]), and allele-specific super-suppressor mutants have been 
reported for the yeast Saccharomyces cerevisiae ([72]). We performed a 
PCR-based pilot screen for Wolbachia candidate tRNA suppressor 
mutants that would be able to rescue both ochre mutations in 
wCerl but did not find any (data not shown). It can hence not yet 
be explained why unexpectedly high frequencies of nonsense 
mutations occurred in two Wolbachia housekeeping genes of both 
Wolbachia strains. We propose, however, several ideas that might 
explain the compensation of these mutations. First, the 'codon 
capture model' allows a bacterial codon that has fallen to low 
frequencies to be reassigned without triggering fitness implications 
([73,74]). In the case of the ochre mutation in wCer2, TGA would 
be re-coded into a synonym and hence not affect protein length. 
Alternatively, in concert with our theory of co-existing roCer2 
haplotypes in the population of one single host, it is possible that 
one haplotype carries the nonsense mutation whereas another one 
does not. A fully functional wCer2 haplotype could then 
potentially compensate the mutation in the non-functional 
haplotype. Generally, bacteria are assumed to be monoploid i.e., 
they carry only one copy of a circular chromosome. Recent 
publications have demonstrated that this is not necessarily case. 
[75] have shown that Neisseira gonorrhoeae are polyploidy and 
carry three genome copies in average. [76] have added striking 
new findings by stating that monoploidy is not typical for bacteria. 
In contrast, polyploidy is very common in proteobacteria with up 
to even 20 genome copies. Following these interesting findings, it 
might be possible that Wolbachia also contain more than one 
genome copy per cell. If those copies are different, i.e., one carries 
the mutation and the other does not, the Wolbachia sequence 
polymorphism we detected in this study can be explained. 
Alternatively, but highly unlikely, the formation of paralogues 
via intrachromosomal duplications of the three Wolbachia genes 
coxA, JtsZ and gatB can be employed. Finally, it could also be 
speculated that alternative genetic codes might support compen- 
sating the nonsense mutations. Contrasting to the bacterial code, 
the TAA triplet does for example not lead to a termination signal 
in the ciliate, dasycladacean and Hexamita nuclear genomes as 
well as in the alternative flatworm mitochondrial code. The TGA 
codon can be compensated by even nine alternative codes (source: 
www.ncbi.nlm.nih.go). So far, we do not know how stably these 
mutant haplotypes are maintained within the wCec population but 
an ongoing deep sequencing project of the «/Cer2 genome will 
allow us to screen, in detail, for the presence and maintenance of 
these mutations. 

4.4 Maintenance and frequency shifts of diversity after 
trans-infection 

Multiple Wolbachia strains can coexist within single host 
individuals e.g. ants ([77]) and tephritid fruitflies ([78,39]). 
However, simultaneous persistence of more than one Wolbachia 
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strain within a single host raises the question as to whether inter- 
strain competition for survival and stable persistence does occur. 
In the case of the observed strain replacement in RC20, inter- 
strain competition between wCerl and wCer2 is expected. We 
collected evidence that the initially common wCer2 was no longer 
traceable in G168 after the trans-infection event, supporting the 
idea of strain replacement following inter-strain competition. 
Recent studies demonstrated that Wolbachia infections can occur at 
extreme low titer levels. The persistence of such natural low titer 
Wolbachia infections have been reported in bark beedes ([33]), 
neotropical Drosophila species ([17]), aphids ([79]), D. simulans 
([80]) and tsetse flies ([81]). We hence assume that »Cer2 density is 
extremely variable, making detection of the symbiont impossible 
even by using highly sensitive techniques. However, the artificial 
double infection in the new host background might have triggered 
ongoing inter-strain competition. An explanation for iwCerl 
displacing prevalent a>Cer2 in RC20 might be a negative 
symbiont-host productivity correlation (see 4.5). Beaumont and 
colleagues recently reported on the experimental evolution of bet 
hedging strategies in bacterial populations ([82]). Bet hedging is 
defined as stochastic switching between phenotypic stages ([83,84]) 
in order to facilitate persistence in fluctuating environmental 
conditions. The new host background, representing changed 
environmental conditions for the symbionts, could have led to the 
evolution of a bet hedging-like strategy in RC20, switching 
between two Wolbachia variants. Switch from »Cer2 to oderl is 
correlated with enhanced fecundity in RC20, and it cannot be 
ruled out that such adaptive bet hedging results in switching back 
to o/Cer2 as main infection variant. 

In RC33 we found clear co-existence of wCerl and a>Cer2, 
obviously not subjected to inter-strain competition. Densities of 
both strains seem to be equal in this system, suggesting 
competition-free co-existence. This was already shown for the 
parasitic wasp Leptopilina heterotoma where density of different 
Wolbachia strains was not affected by the presence of other strains 
([85]). Stable multiple infections with Wolbachia were reported 
from D. simulans ([86]), R. cerasi ([39]), and Ae. albopictus hosts ([87]). 
The idea of stable co-existence without inter-strain competition of 
the bacteria is supported by significantly increased rates of 
maternal transmission determined in this line in comparison to 
the rate evaluated shortly after microinjection. This is in 
agreement with our study where wCer2 in D. simulans exhibited 
an initial prevalence of 65% ([29]), and has now reached almost 
complete transmission (95%; data not shown). 

4.5 Trans-infection events claim reproduction costs in 
novel hosts 

[88] demonstrated that aMelPop triggers severe phenotypic 
changes such as decrease of fecundity in the mosquito Ae. albopictus. 
The authors report a clear correlation between host phenotype 
and the endosymbiont Wolbachia. Similar to the situation in 
mosquitoes, we observed that fecundity of trans-infected D. 
simulans lines is affected by the artificially introduced Wolbachia 
infection when measured ten generations after microinjection 
([29]). Data obtained from an ovary screening assay more than 
150 generations post microinjection suggested that female flies 
were still not adapted to the infection. Compared to the uninfected 
D. simulans Noumea TC strain, 80% of the trans-infected RC lines 
displayed decreased ovary sizes. Most interestingly, ovaries of 
RC20 females were significantly larger, as they contained more 
mature eggs than the other RC lines and slighdy larger than the 
uninfected control. The fecundity of RC20 seemed to change from 
very poor at the beginning to enhanced in later host generations in 
the course of this study (DS, personal observation). Hence we 



might observe a correlation between female fecundity and 
Wolbachia infection in RC20. We have tracked a switch in 
o>Cer2/o)Cerl prevalence in RC20 that most likely occurred 
between generations 150 and 167. It might be possible that this 
change of female fecundity is correlated with a Wolbachia strain 
switch in this line. Although we have no direct evidence yet, we 
speculate that a>Cer2 might be negatively correlated with female 
fecundity in this special case of RC20. This line was reported as 
mono-infected with roCer2 in earlier passages ([29]), later 
diagnosed as wCerl&2 double-infected (this study), and eventually 
the aCerl infection has outcompeted oCer2, since the latter strain 
was no longer traceable by PCR in later generations (Figure IB). 
In order to determine when this shift from wCer2 to twCerl took 
place, we analyzed RC20 DNA extracts from seven, randomly 
picked, non-consecutive generations between G140 (beginning of 
this study) and G168 (see time course in Figure 1A). We found a 
switch from z«Cer2 to <»Cerl during a transition period between 
generations F150 and F167, followed by replacement of wCer2 by 
i*Cerl in F168. However, to direcdy prove the correlation 
between wCer2 and host fecundity, and to rule out that the 
microinjection-caused bottieneck did not just lead to accumulation 
of negative effects, further experiments are needed. Re-evaluation 
of our data in a homogenized host nuclear background obtained 
through outcrossing the RC20 line will allow for better analysis of 
OjCer-triggered fitness costs in the host. 

Conclusion 

In this study we aimed at testing if artificial symbiont transfer 
triggers structural rearrangements, and acceleration of SNPs in the 
symbiont genome. Analysis of mobile genetic elements within 
Wolbachia did not reveal rearrangements after arrival of the 
symbiont in the recipient hosts. By assessing SNP frequency in 
three essential Wolbachia genes before and after microinjection, we 
determined that the purifying constraint operating on these loci is 
hardly relaxed after more than 150 host generations. Instead of 
tracing new mutations upon transfer in the recipients, we 
discovered ancestral polymorphic sites in symbiont genes deriving 
from the donor, pinpointing that both wCerl and <x>Cer2 exhibit 
ancestral and cryptic sequence polymorphism in its original host R. 
cerasi. We further uncovered multiple strains in D. simulans lines 
that were previously typed as singly infected. This may have been 
due the co-existence of Wolbachia strains, where one of these 
persisted at low titer and thus had escaped standard detection 
techniques. We demonstrated that infections by multiple strains 
are prone to shifts in strain prevalence upon artificial host transfer. 
This reflects the population-like structure of Wolbachia within and 
between different hosts and thus will have consequences for 
symbiont population dynamics. Persistence of cryptic multiple 
infections after transfer from a multiply infected donor, captures 
the importance of studying, in detail, the integrity of Wolbachia 
infections prior to application as tools in modern pest and disease 
control management. 

Supporting Information 

Data SI Extended methodology for RFLP-mapping via 
genomic Southern blot analysis. Detailed information about 
RFLP mapping can be found in Data SI. 
(DOCX) 

Figure SI Wsp-PCR. of antibiotic treated RC20, -50 and 

R. cerasi. Upper lane: Generation Fl from 0.03% tetracycline- 
treated RC20 (1-2) and RC50 (3-4); generation F2 in same order 
(5-8). RC 50 harbors Wolbachia at higher densities than RC20 
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resulting in still bright signals upon treatment in Fl and highly 
significant reduction in F2. Positive control (Wolbachia-miected 
DSR) is shown in 9. Lower lane: untreated R. cerasi larva (1); 
antibiotic-treated R. cerasi larvae from 0.02%, 0.025%, 0.05% 
tetracycline (2-4); antibiotic-treated R. cerasi adults from 0.02%, 
0.025% tetracycline (5-6); untreated R. cerasi adult (7). Negative 
controls [Wolbachia-urimitcttd D. simulans Noumea TC) are shown 
in 8 and 9. 
(TIF) 

Figure S2 RFLP analysis of wCer2 infection in single 
trans-infected Drosophila (RC) via Southern hybridiza- 
tion. (A) Hybridzation with ZSJ-probe (Insertion Sequence 
Element), on //mdlll-digested total DNA of D. melanogaster 
Harwich (aiMel; WoP os ) as reference and trans-infected RC lines 
RC33 (1), RG45 (2-3), RG50 (4-5), RC21 (6). Black asterisks 
indicate two fragments that are not present in the aMel reference 
and are thus considered diagnostic for wCer2. i/mdIII-lambda 
bands 9400, 6600, 4300, 2300, 2000. (B) Number of fragments in 
characteristic RFLP-fingerprints of aJVIel and wCer2 using six 
different Wolbachia marker. "Expected numbers of fragments after 
in silico analysis of the annotated wMel genome (indicated by 
asterisk; GenBank accession number NC_002978.6), and b exper- 
imentally determined fragment numbers obtained from RFLP 
analysis. Number of fixed fragments refers to fragments that are 
present in the reference «Mel and not considered diagnostic for 
»Cer2. (C) Hybridization with IS3 (D) ISNew (E) FJV7S-141 (F) 
upper panel: K7V7S-144, lower panel: wsp. Abbreviations: M 
Marker (Lambda-DNA digested with HindllT). 
(TIF) 
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